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Abstract 

We study spherical, charged and self-similar distributions of matter 
in the diffusion approximation. We propose a simple, dynamic but phys- 
ically meaningful solution. For such a solution we obtain a model in 
which the distribution becomes static and changes to dust. The collapse 
is halted with damped mass oscillations about the absolute value of the 
total charge. 

1 Introduction 

Several authors have considered charged distributions of matter El ^ |S1 

IS »L ^ JiL lUj I although it is well known that astrophysical objects are not 
significantly charged. Nevertheless, in some stages of the gravitational collapse 
even a small amount of charge can change the final state of the body. Some in- 
teresting features of charged collapsing matter justify any effort to get physical 
insight by studying this problem. For instance, naked singularities can be pre- 
vented ^21; the final geometrical structure left over after the complete collapse 
of a spherically symmetric charged source and of a chargeless rotating star are 
similar Cauchy horizons, gravitational repulsion and perhaps traversable 
wormholes are also possible 

If the mathematical treatment is simplified, the evolution of a charged dis- 
tribution of matter can be followed by the Einstein-Maxwell equations. In this 
paper we explore the self-similar gravitational collapse of charged spheres in the 
diffusion approximation. It is well known that dissipation due to the emission 
of massless particles is a characteristic process in the evolution of massive stars. 
The only plausible mechanism to carry away the bulk of binding energy of the 
collapsing star, leading to a neutron star or black hole, is neutrino emission 
[Tl) . It seems clear that the free-streaming process is associated with the initial 
stages of the collapse, while the diffusion approximation becomes valid toward 
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the final stages. The junction conditions at the boundary of a charged and 
dissipative sphere have been considered with outgoing heat flow and radiation 
flux [151 On the other hand, the field equations admit homothetic mo- 
tion |17lll8[[Tmi2(J) . Applications of homothetic similarity range from modeling 
black holes to producing counterexamples to the cosmic censorship conjecture 
inilSSEniElOSlEniOZlllHl- in particular, homothetic charged and isotropic 
(or anisotropic) fiuids have been studied pUll5ni f7|. 

We observe in the literature that much work has been done under static 
conditions and dusty charged matter. Also authors make additional assumptions 
such as equations of state or relationships between metric variables jSOl |^ . In 
this paper we obtain a dynamical model from a simple solution to the homothetic 
motion, without any additional supposition except spherical symmetry and self- 
similarity. In section 2 we write the field equations and the junction conditions. 
The equations at the surface of the distribution of matter are presented in section 
3. In section 4 we write the symmetry equations to describe self-similarity. In 
section 5, we show an example from a simple solution which we discuss in the 
last section. 



2 Field equations and matching 

Let us consider a non-static distribution of matter which is spherically symmet- 
ric and consists of charged fluid of energy density p, pressure p, electric charge 
density a and radiation energy flux q diffusing in the radial direction, as mea- 
sured by a local Minkowskian observer comoving with the fluid. In radiation 
coordinates [321 metric takes the form 

ds"^ = e^'^ (^du^ + 2du dr^ - [dO^ + sin 6^ d(f)^^ , (1) 

where /3 and V are functions of u and r. Here u is the time-like coordinate: 
in the flat space-time u is just the retarded time, so that surfaces of constant 
u represent null cones open to the future; r is a null coordinate (g^j. = 0) such 
that surfaces of constant r and u are spheres; 9 and (f> are the usual angular 
coordinates. Also, in this paper we use relativistic units {G = c = 1). 

The Einstein field equations, G — — SttT, are considered with the energy- 
momentum tensor T = {p + p)v^\-—pg+q^v+v^q+E, where v and q are 
the four-velocity and the heat flux four-vector, respectively, which must be 
orthogonal; E is the electromagnetic field energy-momentum tensor constructed 
with the Maxwell field tensor F as is usually done. The Maxwell field equations, 
d*F= 47r*J and dF= 0, are coupled minimally with gravitation, where J= av 
is the electric current four-vector. 

Because of the spherical symmetry, only the radial electric field F"'' — —F™ 
is non- vanishing. If we define the function C{u,r) by the relation 

i?^«- = Ce-2/3/r2, (2) 
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the inhomogcncous Maxwell equations become 

C,r = 47rrV"e2^ (3) 

and 

C„ = -inr^re^f", (4) 

where the comma subscript represents partial derivative with respect to the 
indicated coordinate. The function C{u, r) is naturally interpreted as the charge 
within the radius r at time u. Thus, the inhomogeneous field equations entail 
the conservation of charge inside a sphere comoving with the fluid, expressed by 

v(C) = 0. (5) 

Let uj be the velocity of matter as seen by a Minkowskian observer moving 
at —Lu with respect to the local comoving frame; the matter velocity in radiation 
coordinates is then given by 



dr V UJ 



(6) 



du r 1 — UJ 
Introducing the mass function by 

V = e^f^ir ~ 2m{u, r) + C^/r), (7) 

we can write the Einstein equations as 

p + puj^ 2ujq e^^^{CC,u/i — rh^u) fh^r — CC^r/r 



1-^2 l-uj^ A'!rr{r ~2m + C^/r) Anr"^ 



(8) 



p — pUJ f 1 — UJ^ rhr — CC^r/l^ 



l+uj \l+uj J Anr^ 
^ " ""Vp + P) - 2 f ^) q=^{r-2m + C'/r), (10) 



1 + UJ J \1 + uj J 2'Kr 

P = -^P,ur e-^^' + ^(1 - 2m/r + C"" /r''){2(3^„ + ^I3\ - f3,r/r) 
An 8tt 

+ [3/3 r(l - 2m r) -mrr] + ^(2CC^/r - C^/r^) 

+ F^(^r + CC\rr - 2CC.r/r). (11) 

We describe the exterior space-time by the Reissner-Nordstrom-Vaidya 
metric 

dsl =(^1- + ^^du^ + 2du dr - {dO^ + sin ^9 dcjj^) , (12) 

where m{u) is the total mass and Ct the total charge. It can be shown that 
the junction conditions to match the metrics Q and p2|l . across the moving 
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boundary surface r = a(it), are equivalent to the continuity of the functions 
m{u,r) and I3(u,r) across the boundary, that is, •m{u,a) = •m[u) and P{u,a) = 
0, and to the equation 



- ^['-- + #)^'™ - ^ + i/" = 0' (1^) 

where the subscript a indicates that the quantity is evaluated at the surface 
r = a{u). We have used the continuity of the radial electric field through the 
boundary assuming no surface free charge density, yielding C{u,a) = Ct- It 
should be mentioned that the discontinuity of the pressure at the boundary, 
?5a = 9a, is a direct consequence of (jHl, ®, HlOU and the junction conditions 

123113 HSI 



3 The surface equations 

One of the surface equations is just © evaluated at r = a, which takes the form 

a^(l-'- + ^)^, (14) 

where {TJ has been used as well as the junction conditions for m, (3 and C . It 
is convenient to scale the variables by the initial mass r7i(0), such that 

/I - " A-f - " '^^ n 
A — — -— ; M = — -— ; — -— — > u\ — > Gt- 

to(0)' m(0) m(0) ' m(0) 

Defining the surface potential F and U, by 

and 

n ^ (16) 

1 - Wa 



p4|l can be written in the form 

i = (17) 

which is the first surface equation. 

The second surface equation relates the total mass loss rate to the energy 
flux through the boundary surface. It has been shown \6'6\ that this equation 
can be written as 

. 2QF+(l-f)i CU 

P= J (18) 

where 

Q^i2n-l){ATrr^q)r=a. (19) 
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The third surface equation is the charge conservation law given by (O , evaluated 
at the boundary 

The fourth surface equation is model dependent. This correspond to the 
Bianchi identity T/*^ = 0, given by 

dp P + P ( , -20 ( P + P \ 

dr {l-2m/r + Cyr^)\''''^ r^J \1 - 2m/r + C^r^ J 

--{P-P) + ^ (21) 

where p = (j> — ^p) /(I +w) — (1 — ui)q/{l -\- uj) and p = {p — top) /(I + uj) — 
(1 — (jj)q/{l + Lo). Equation (|21|l is the generalization of that of Tolman- 
Oppenheimer-Volkov for non-static and charged radiative situations j23| (see 
Reference 0J for the static case). 



4 Self— similar space— time 

Self-similarity is invariably defined by the existence of a homothetic Killing 
vector field J^. We shall assume that the spherical distribution admits a one- 
parameter group of homothetic motions. A homothetic vector field on the man- 
ifold is one that satisfies £(^g =2ng on a local chart, where n is a constant on 
the manifold and £ denotes the Lie derivative operator. If n 7^ we have a 
proper homothetic vector field and it can always be scaled so as to have n = 1; 
if n = then ^ is a Killing vector on the manifold |34[ I35| . So, for a constant 
rescaling, ^ satisfies 

£iS = 2g (22) 

and has the form 

^^Aiu,r)du + X{u,r)dr. (23) 
The homothetic equations reduce to 

e(X) = Z^{Z), (24) 

ay) = 0, (25) 

where A = r, A = A(u), X = m/r, Y = Ae^^/r and Z = C/r. Therefore, 
X — A(C), Y = Y{C) and Z = Z{() are solutions if the self-similar variable is 
defined as 

CEEre-''*'/^. (26) 

In the next section we propose a simple solution which evolves toward staticity, 
rendering an inhomogeneous and dusty fluid sphere. 
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5 A simple example 



The simple solution rh — mr/a, e^^ — r/a and C — CttIo, satisfies the ad- 
ditional symmetry equations (|24() and the charge conservation equation 
H2U|I . the continuity of the radial electric field and the continuity of the first 
fundamental form. Now, with these solutions and the junction condition H13|l 
(continuity of the second fundamental form), f2 is determined by 

(.7, 

The heat flow at the surface is obtained from the conservation equation H21|l 
evaluated at the surface, resulting in 

Q=^-^^l^{\-2F-CllA\ (28) 

This last equation, together with (|27ll . must be taken into account to integrate 
numerically the equations p7|) and (|18|) . Using a standard Runge-Kutta algo- 
rithm and the initial conditions A(0) — 3.50 and m(0) = 1.00, we study the 
effect of charge on collapse. Once the boundary evolution and its energetic are 
determined, we calculate the physical variables from the field equations. Figures 
1-7 show the results obtained for our simple solution. We shall discuss them in 
the next section. 



6 Discussion 

Figure 1 displays the evolution of the surface for different values of the total 
charge. It is clear that the increase of the total charge favors the collapse in a first 
stage of the evolution. Later the collapse is halted, with damped oscillations 
rendering a distribution which is less compact for greater total charge. We 
show in Figure 2 how the Bondi mass decreases and oscillates about the total 
charge (given positive) until reaching the same value of Ct- Figures 3-6 sketch 
the pressure p, the density p, the heat flow q and the matter velocity dr/du. 
Observe that the whole sphere of fluid becomes dust (p = at all points) and 
inhomogeneous when staticity is reached. It is interesting to note that the 
cooling proceeds with emission and absorption of energy and consequently all 
the shells bounce and contract as a unit, over and over, until reaching staticity 
in the whole sphere. Also it is striking how the ratio p/p is a function only of 
the Bondi time, i.e., it is the same at any point of the charged distribution (see 
Figure 7). 

We would like to stress that our solution, in spite of its simplicity, behaves 
very well when one takes into account the results reported by other authors. 
The sphere collapses and rebounds as a unit jl3| . over and over, until reaching 
staticity. Equilibrium configurations are possible, with no necessity of inter- 
nal pressure; moreover, electric charge halts the gravitational collapse 0|. 
The spheres of charged matter can oscillate and the final static configuration is 
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reached when the total mass is equal to the total charge |2I, that is, the charged 
spheres of dust in equilibrium belong to the interior Papapetrou-Majumdar 
class ini- 

Finally, we make special mention of the Herrera and Ponce de Leon ho- 
mothetic models [7]. For a null radial pressure at the surface (as a boundary 
condition) they find a charged dust. Otherwise, the distribution is infinitely 
extended. From an altogether physical point of view, perhaps it is possible to 
find finite perfect (and self-similar) fluid sources with zero pressure at the sur- 
face, considering free-streaming as the inner transport mechamism. In general, 
Herrera and Ponce de Leon find that the total charge is less than the total mass. 
The equality will only be valid in the case of homothetic charged dust spheres. 
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Figure 1: Radius a as a function of the time u for A{0) = 3.5, to(0) = 1 and 
for different values of the total charge Ct- 0.01 (dot-dashed line); 0.05 (dashed 
line); 0.09 (dotted line); 0.2 (solid line). 



9 




Figure 2: Bondi mass at the surface m as a function of time u for A{0) = 3.5, 
m(0) = 1 and for different values of the total charge Ct- 0.01 (dot-dashed 
line); 0.05 (dashed hue); 0.09 (dotted hue); 0.2 (solid hue). Observe that the 
final mass in each case is equal to the total charge. 
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Figure 3: Pressure p as a function of time u for ^(0) = 3.5, to(0) = 1 and 
Ct = 0.09 at different points: r/a = 0.25 (dot-dashed line); r/a = 0.33 (dashed 
line); r/a = 0.50 (dotted line); r/a=1.00 (solid line). 
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Figure 4: Density p as a function of time u for A(0) = 3.5, m(0) = 1 and 
Ct = 0.09 at different points: rja = 0.25 (dot-dashed line); r/a = 0.33 (dashed 
line); rja = 0.50 (dotted line); r/a=1.00 (solid line). 
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Figure 5: Heat flow q as a function of time u for A{0) = 3.5, m(0) = 1 and 
Ct = 0.09 at different points: r/a = 0.25 (dot-dashed line); r/a = 0.33 (dashed 
Hne); r/a = 0.50 (dotted line); r/a=1.00 (solid line). 
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Figure 6: Matter velocity dr/du as a function of time u for ^(0) = 3.5, m(0) = 1 
and Ct = 0.09 at different points: r/a — 0.25 (dot-dashed line); rja — 0.33 
(dashed line); rja = 0.50 (dotted line); r/a=1.00 (solid line). 
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